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Electromagnetic Form Factors and Charge Densities From Hadrons to Nuclei 



Gerald A. Miller 
Department of Physics, University of Washington 
Seattle, Washington 98195-1560 

A simple exact covariant model in which a scalar particle ^ is modeled as a bound state of two 
diflterent particles is used to elucidate relativistic aspects of electromagnetic form factors F{Q^). 
The model form factor is computed using an exact covariant calculation of the lowest-order tri- 
angle diagram. The light-front technique of integrating over the minus-component of the virtual 
momentum gives the same result and is the same as the one obtained originally by Gunion et al. by 
using time-ordered perturbation theory in the infinite-momentum-frame. The meaning of the trans- 
verse density p{b) is explained by providing a general derivation, using three spatial-coordinates, 
' of its relationship with the form factor. This allows us to identify a mean-square transverse size 

^— I (6^) — J d?h p{h) = —'^-^^{Q^ = 0). The quantity (6^) is a true measure of hadronic size because 

of its direct relationship with the transverse density. We show that the rest-frame charge distri- 
{^JQ' bution is generally not observable by studying the explicit failure to uphold current conservation. 

^ , Neutral systems of two charged constituents are shown to obey the conventional lore that the heavier 

one is generally closer to the transverse origin than the lighter one. It is argued that the negative 
central charge density of the neutron arises, in pion-cloud models, from pious of high longitudinal 
momentum that reside at the center. The non-relativistic limit is defined precisely, and the ratio of 
the binding energy B to the mass M of the lightest constituent is shown to govern the infiuence of 
relativistic effects. We show that the exact relativistic formula for F{Q^) is the same as the familiar 
one of the three-dimensional Fourier transform of a square of a wave function for very small values 
of B/M, but this only occurs values of B/M less than about 0.001. For masses that mimic the 
quark-di-quark model of the nucleon we find that there are substantial relativistic corrections to the 
O ' form factor for any value of . A schematic model of the lowest s-states of nuclei is developed. 

^ , Relativistic effects are found to decrease the form factor for light nuclei but to increase the form 

C . factor for heavy nuclei. Furthermore, these lowest s-states are likely to be strongly influenced by 

relativistic effects that are order 15-20%. 
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I. INTRODUCTION 



' The text-book interpretation of nucleon electromagnetic form factors is that their three-dimensional Fourier trans- 
^\ \ forms are measurements of the charge and magnetization densities. This interpretation is deeply buried in our thinking 
■ and continues to guide intuition as it has since the days of the Nobel prize- winning work of Hofstadter[l| . Nevertheless, 
the relativistic motion of the constituents of the system causes the text-book interpretation to be incorrect. 

The preceding statement leads to a number of questions, the first being: Is the statement correct? If correct, how 
relativistic does the motion of the constituents have to be? Why is it that the relativistic motion of the constituents 
and not that of the entire system that causes the non-relativistic approach to fail? It is probably true that the answers 
to these questions are displayed within the existing literature. However, obtaining general clear answers has been 
sufficiently difficult that posing even the first question of this paragraph would not lead to a unanimous answer by all 
professionals in the field. 

This paper offers the strategy of using a simple model, a generalization of the cjT' model used by Weinberg 0] to 
illustrate advantages of using the infinite momentum by choosing frame that was used by Gunion et al. ^] to explore 
form factors and hadronic interactions at high-momentum transfer. We take the interaction Lagrangian density to 
be of the form g'^4)£, in which all of the fields are bosons. The 4* particle of mass M represents the bound state of the 
two different constituents 0,^ of masses TOi,m2 respectively. Thus the represents the hadron or nucleus with the 
0,^ representing the quark or nucleonic constituents. Mass renormalization effects are ignored. We can choose either 
or both of the constituents to be charged and thus discuss charged and neutral ^' particles. 

The motivation to pose questions regarding the meaning of electromagnetic form factors at this point in time arises 
from recent experimental work, especially the discovery that the ratio of the proton's electric to magnetic Sachs form 
factors Ge/Gm drops rapidly (please see the reviews and from our recent finding Q, based on measurements and 
the use of the transverse density that the charge density at the neutron's center is negative. The nucleon transverse 
density p{h), the two-dimensional Fourier transform of -Fi is the infinite momentum frame charge density [6] located at a 
transverse separation b from the center of transverse momentum^, d, fiot . This quantity has a direct relationship to 
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FIG. 1: Feynman diagram for the form factor with the photon coupling to the <j) particle of mass mi. The initial and final 
hadron $ carry momentum, P and P + q. The ^ is a spectator. 

matrix element of a density operator. The usual three-dimensional Fourier transforms of and Gm do not because 
the initial- and final-state nucleons have different momentum, and therefore different wave functions. This is because 
the relativistic boost operator that transforms a nucleon at rest into a moving one changes the wave function in a 
manner that depends on the momentum of the nucleon. However, we expect that there are non-relativistic conditions 
for which the text-book interpretation is correct. We aim to explore those conditions by choosing appropriate values 
of the masses mi, TO2, M . 

Here is the outline we follow. The form factor for the situation in which the 5* and carry a single unit of charge, 
but the ^ is neutral, is computed using an exact covariant calculation of the lowest-order triangle diagram in Sect. H. 
This is followed by a another derivation using the light-front technique of integrating over the minus-component of the 
virtual momentum in Sect. HI that obtains the same form factor. This is also the result obtained originally by Q by 
using time-ordered perturbation theory in the infinite-momentum-frame IMF. Thus three different approaches yield 
the same exact result for this model problem. Any approximation that does not yield the same form factor is simply 
not correct. The asymptotic limit of very high momentum transfer is also studied. The next section (IV) explains 
the transverse density of the model, its central value, a general derivation of its relationship with the form factor 
using three dimensional spatial coordinates and the meaning of hadronic radii. Section V displays the spatial wave 
function in terms of three spatial coordinates. Section VI shows that the rest-frame charge distribution is generally 
not observable. Section VII is concerned with the question of whether neutral systems of two constituents obey the 
conventional lore that the heavier one is generally closer to the origin than the lighter one. The non-relativistic limit 
is defined and applied to a variety of examples in Section VIII. The exact formula for the form factor morphs into the 
familiar one of the three-dimensional Fourier transform for sufficiently large values of the constituent mass divided 
by the binding energy of the system. Examples that are motivated by the pion, deuterium, nucleon and heavy nuclei 
are provided. This work is summarized in Section IX. 



The model Lagrangian density is given by g^^cj) ^ where and f represent three scalar fields of masses M,mi 
and m2 respectively and g is a coupling constant. One can take two or three of these fields to carry charge to make 
up a system of definite charge (including the case when the hadron is neutral) . The effects of mass renormalization 
are not considered here because we aim to use a simple model to provide easily calculable examples and illustrate 
specific points. The condition mi -I- m2 > M is used to insure that the hadron ^E" is stable. 

We start with the situation in which the '5, (f> each carry a single positive charge and ^ is neutral. The form factor 
F{q'^) for a space-like incident photon of four-momentum qf^ [q^ < 0,Q^ = —q^), incident on a target 4* of four- 
momentum interacting with the of mass mi is given, to lowest order in g, by the single Feynman diagram of 



II. EXACT FORM FACTORS USING A SIMPLE MODEL 
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Fig. 1 . We take the model electromagnetic current (in units of the proton charge) as given by 



and find 



(27r)4 [k^-mi + ie] 



{2k^ + q") 



{{k + qY - mj + ie) {{P - k)^ - + ie) ' 



(1) 

(2) 
(3) 



Proceed by combining the denominators using the Feynman procedure and shifting the origin of the convergent 
integral to find 



dx 



(2^)^ 7o " Jo 

where = Q^y{y + x — 1) + APx{l — x) — m\{l — x) — m\. Use 

d^/t 1 



dy 



q^'il - 2y) + 2Pf'x 
[K^-M^+ie]3 ' 



TT^ 1 1 



(27r)4 (k2 -M^ + ie)3 



to find 
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dx 



dy- 



(27r)4 2 7W2 
q''{l-2y) + 2Pfx 



^y{y + a; — 1) + M^x{l — x) — ml{l — x) — m\x ' 



(4) 



(5) 



(6) 



The integral over y can be done in closed form with the result 

:rTanh"^ 



F{Q'^){2P^ + q^ 
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■(q^ + 2P^) / dx 



^Ax m|+4(l-a;)mf-x(l-2;)M2 + (l-a;)2 



/(^.jAx ml + 4(1 - x)ml - 4a;(l - x)NP + (1 - a;)2Q2 



(7) 



The above expression shows that current conservation is satisfied and that the form factor can be obtained from any 
component of the current operator . The final result for the form factor is 



xTanh 



dx 



yjAx ml+A(\-x)m\~x{l-x)M'^ + {\~xYQ'^ 



47^2 7o ^Jl^^^x ml + 4(1 - x)ml - 4x(l - x)AP + (1 - a;)2g2 ' 



(8) 



This closed form expression is the key result of this paper. It can be modified to describe a variety of different physical 
situations. 



III. INFINITE MOMENTUM FRAME/LIGHT FRONT REPRESENTATION 

We derive the light front representation by starting with the form factor of Eq. ([3]) and integrating over k^ . This 
procedure is simplified by choosing g+ = 0, so that = @| and evaluating the + component of the electromagnetic 
current operator. In the present Section the convention is that — A'^ ± A'^ for the four- vector A^. Then Eq. Q 
becomes 



2k+ 1 1 



(fc2 - m^l + ie) {(k + qf - mj + ie) ((P - fc)2 - AP + ie) 



(9) 



= _ ■ 2 f d^k 2k+ 1 1 1 

J (27r)4 fc+2(p+ _ _ ^ ^) _ ((k+q)^+W) ^ ^) (p- _ k- (£^^1!+^ + ^ 

If we integrate over the upper half of the complex k^ plane we find a non-zero contribution only for the case 
Q < k^ < P+ . Carrying out the integral leads to 

')P+F(n^\ = [ rfi\c f I ], (-1-1) 

) (27r)3 / " / k+(P+-k+)p- k^+mf (P-k)f+m| p_ (k+q)^+mf (P-kj^+mj 
^ ^ ^ V J I- ^ p+-k+ ^ feT P+-k+ 
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Next change variables by defining 



so that 



Fm = 
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dx 



1 



X(l - a;) p+p- _ k2+rn^ _ (P-k)^+m| p+p_ _ (k+q)^+m^ _ (P-k)2+mi 

X 1 — a; a; 1— a: 



2(27r)3 

Further define the relative transverse momentum 

K = (1 - x)k - x(P - k) = k - xP, 
so that the form factor can be re-expressed as 



2(27r)3 



dx 



1 



X{1 - x) f^j2 _ ^"+"^1 _ H'^+ynj j^j2 _ (^+(l-^)q)"+"^i _ (/t+(l-x)q)^+ri 



(12) 



(13) 



(14) 



(15) 



This is the expression obtained in Ref. by using time-order-perturbation theory in the infinite momentum frame. 
Integration over k~ leads to the same result for this example. 

It is useful to re-express the result Eq. ()15|) in terms of a wave function ip with 



ij{x,K) = g[M^ 



'21-1 



l-X 



In that case 



2(27r) 



d'^n 



dx 



-ij*{x, K,+ {I - x)q)ip{x, k), 



(16) 



(17) 



/o a;(l - x) 
as found in Ref. 

The integration over /t is convergent so we carry out the integration over k by combining the propagators and 
shifting the origin. This gives 



F(Q2 
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2(27r)3 



dx / dz 



x{l — x) 



(1 — x)mi + xm\ — x{l — x)K'P + (1 — xYa^z{l — z) 



(18) 



The integral over z can be done with the result that 



PiQ') -i-, dx 



xTanh 



^yix m|+4(l-a;)mf-2;(l-a;)M2 + (l-a;)2Q2 



47^2 Jo ^/O^^ix ml + 4(1 - x)ml - 4x{l - x)M^ + (1 - x)'^Q'^ 



(19) 



This is the same as our previous exactly computed result, Eq. ([8]). Thus evaluation in the infinite momentum frame 
or the equivalent (for this model) light front technique of integration over k~ yields the exact result. 



A. Asymptotic Behavior of the Form Factor 

The limit of very high is of considerable interest. One wants to see how the quark counting rules emerge from 
an exact calculation, even if the model is very simple. To this end we note, that the integral Eq. (|19|) can be evaluated 
exactly in the limit that toi = m2 = m with AI — 0. Then measuring Q = \J in units of m {Q/ni — > Q) we find 



so that 



^ ^2 log2(i(g(v/Q^+Qj +2jj +8 v/Q^log(i(Q(v/Q^+Qj+2j 
^^^^^4^ 8Q5 2Q' ~ 

ilog2(i)-Hlog(^)+l log(^)-l 
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Thus the leading asymptotic behavior is 



hm F{Q') ^ ^i^. (22) 



Thus the power-law fall-ofF expected from the quark-counting rules appears, but it is modified by the presence of 
the logarithms. This behavior is not associated with taking M to zero because in all cases we have mi -I- m2 > M 
as required for the particle to be stable. Thus the asymptotic behavior {Q^ ^ rnl^rn^) associated with Eq. (|22p is 
expected to be universal for this model. Note however, from Eq. (j21[) that the approach to this asymptotic form is 
very slow. 



IV. ELECTROMAGNETIC FORM FACTORS MEASURE TRANSVERSE DENSITIES AND 

TRANSVERSE RADII 

The expression Eq. (jl7p is noteworthy because the form factor is expressed as a three-dimensional integration 
that involves momentum-space wave functions evaluated at different initial and final momenta. If the {1 — x) factor 
multiplying q were replaced by a constant Eq. (|17p would be similar to the usual expression for the form factor. We 
clarify this comparison by expressing the wave function of Eq. (|16|) and Eq. (|17|) in transverse position space, with B 
canonically conjugate to the transverse momentum variable k: 

^{x, B) = ^=L= / ^e»'^-BV(x, /.), (23) 
^x{l-x) J (27r)^ 

_ \/x{l - x) 
27r 



gKQ{ Jml{l - x) + mix - M'^x{l - x) B), (24) 



with the phase space factor , incorporated in the wave function. Then the form factor Eq. (jlZp can be 

re-expressed as 

^^^'^ ^ 2(hf [ / (25) 

Further simplify by replacing the relative transverse position variable B by the transverse position variable of the 
charged parton bi = b. We have 

B = bi - b2 (26) 
= bi(a;)-f b2(l-x), (27) 
B = b/(l-x), (28) 

where the middle equation sets the transverse center of P+ momentum to zero. Use Eq. (1^51) in Eq. (P5|) to find 



which can be re-written as 



(30) 



with the transverse density p{b) given by 



Pib) - 7- / = tJ;-^ / dx——^K'M mli^ - ^) + m^x - APx{l - x) —-) (31) 



47r7o {l-xY'^' ' l-x" 2(27r)3 7o " (1 - a;) ^ ^ 'l^x 

The transverse density p{h) has been derived previously Q as the integral of the impact parameter generalized 
parton distributions GPD p(x, h) over all values of x. The quantity p{x, b) gives the probability that a quark of 
longitudinal momentum fraction x resides at a transverse position 's', ^ . For the present model 

''^^'^^^ 2(2^)3 ^^l^^ K^i^rnlil-x)+mlx-M^x{l-x) ^). (32) 
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The transverse density is also the integral of the three-dimensional infinite momentum frame density p{x~ , h) over all 
values of the longitudinal position coordinate [ll| . 

The transverse density is directly obtainable from experiment via the inverse Fourier transform of Eq. (|30p provided 
the electromagnetic form factor is measured for sufRciently large values of . The momentum transfer is transverse 
in direction so that information about the longitudinal position or momentum is not available. There is no way to 
use only measured values of F{Q^) to determine p{x, b). 



A. Singular central density 

Before proceeding it is worthwhile to point out that the model central density is singular. This arises as a conse- 
quence of the zero range nature of the coupling. The transverse density p{b) is an integral involving the singular 
function Ko{x) which varies as log i for a; << 1. The question of the singularity of the central transverse density p{b) 
is interesting to the present author because of recent work [ll| showing that for the pion p{b) is likely to approach 
infinity as b approaches zero. We may study the limit as b approaches by using the asymptotic limit of the form 
factor Eq. (|22p . The density for 6 near zero is controlled by the form factor at large values of Q^. We use the inverse 
of Eq. (gni) to write 

limp(6)^ r^'^ln2(g), (33) 

where Qo is a momentum transfer large enough so that so that Eq. (I22p is valid, and e is a fixed positive number 
small enough so that Jq (e) = 1 to any desired precision. Changing variables to u = In Q shows that 

limp(6)^ln3(fe)/3, (34) 
which is the central singular charge density arising from the log^ Q^/Q^ behavior of the asymptotic form factor. 



B. Transverse Charge Density from a more general perspective 

For a spin-0 system, the form factor F{Q'^), Eq. ^ may be computed, in the Drell-Yan DY frame (g+ = 0,(5^ > 
= q^), by using 

The spatial structure of a nucleon can be examined if one uses 0, B [1] ■ The state with transverse center of mass R 
set to is formed by taking a linear superposition of states of transverse momentum: 

\p+,R = 0,)=JV J -^\p+,p,), (36) 

where p, A) are plane wave states and A/" is a normalization factor satisfying \Af\^ J ~ ^- normalization 
of the states is given by 

(p'+,p'b+,p> = 2p+(2^)3j(p'+-p+)<5(2)(p'_p). (37) 

References [H, [l^ use wave packet treatments that avoid states normalized to S functions, but this leads to the same 
results as using Eq. ([36]) . Note however, the relevant range of integration in Eq. ([36]) must be restricted to |p| <C p+ to 
maintain the interpretation of a nucleon moving with well-defined longitudinal momentum [l^ . Thus we use a frame 
with very large p'^ . It is in just such a frame that the interpretation of a nucleon as a set of a large number of partons 
is valid. 

We evaluate the density operator in the infinite momentum frame in which the spatial coordinates are = 
{t — z)/^/2, b, and time is x~^ = {t + z)/V^ = 0. We therefore do not write the x~^ dependence in any function below. 
The infinite momentum frame charge density operator (in units of the proton charge) is given by 

Poo{x-,h) = J+(a;-,b) =0 9+ 0(a;-,b), (38) 
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and the density itself by 



^ ,^ (p+,R = 0|Poo(x-,b) b+,R = 0) 



Use translational invariance in the form Poo (a; ,b) = e*^^^ e 'P '^poo(0)e+'P ''e along with Eg. ([39)1 . Eg. ((36 l) 

and Eq. to determine that 

j dx-poo(x-,b) = ^ y rf2q^^(g2 ^ q2)e-'q b ^ (40) 
Thus one recovers the two-dimensional Fourier transform of Eq. (jSOp . 



C. Mean-squared transverse radii and mean-squared effective radii 



The two-dimensional Fourier transform of Eq. (|30p may be expanded as a power series in Q as 



Jim^F(Q2) = l-^(62), (41) 
where the mean-squared transverse radius (6^) is given in terms of the transverse density as 

(62) ^ / dHb^p{h), (42) 



and a direct relation with the transverse density is evident. In contrast, the usual procedure is to write 



lim F(Q2) = 1 _ ^i?*^ (43) 

Q^— ►0 6 



where we denote i?* the effective mean-squared radius [14] . The quantity i?* has no direct relationship with a density 
unless the system is non-relativistic. Thus we maintain that {b'^) is the basic quantity related to an underlying density. 
However, once the effective mean-squared radius R*^ is determined, the fundamental (6^) is known immediately 
because a comparison of Eq. (|iT|) and Eq. reveals that 



(62) = Ir*\ (44) 



V. WAVE FUNCTION AS A FUNCTION OF THREE POSITION VARIABLES 

The previous section shows that the form factor is simply related to the three-dimensional coordinate-space density 
that depends on [x~ ,h) in the infinite momentum frame. Given the simplicity of our model, we should be able to 
identify a wave function and density. 

The basic idea is that the position variable of a particle is canonically conjugate to the plus-component of the 
momentum. The momentum of the c harg ed constituent is = xP^ , and its canonical longitudinal position variable 
is with [a;~,fc"'"] = i = [x~ ,x\P~^ [15j . The canonical longitudinal position variable for the other particle can be 
taken as —x~. So we can convert the wave function ip{x^ B) of Eq. (|24[) to one expressed entirely in coordinate space. 
We find 

tp{x- ,B) = [ V(a;,B)e"^''^", (45) 

V 27r Jq 

which preserves the normalization condition that F{Q^ = 0) = 1. Note that B is a relative variable and that x~ 
is the variable for the position of the charged constituent. This wave function displays no spherical symmetry-the 
longitudinal x^ and transverse position dependence B are not related. Another point is that the wave function 
Eq. (|45p explicitly depends on the momentum P"*". As P"*" approaches infinity, the value of x~ must be very small to 
prevent x~P~^ from being very large and causing the integral on Eq. (j45|l to vanish. This means that the system can 
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FIG. 2: Two TOPT diagrams for the form factor with the photon coupling to the particle of mass mi. 



be thought of as having a pancake or disc shape. For this reason, the position b = really does correspond to the 
center of the hadron. We shall show below, that in the non-relativistic limit, rotational symmetry emerges. 

The wave function can be computed in closed form for the special case M — 0, mi = m2 — m. Using Eq. (|24p in 
Eq. ([IS]) with the stated parameters leads to the result 



V'(x-,B) 




g Ko(m B) I dx 



\/xil - x) 



cP+a 



2n 



Ko{m B 



(46) 
(47) 



For this simple example, the x^ and B dependence factorizes, showing the explicit violation of rotational symmetry. 
The formula Eq. ((47l) shows also how the spatial extent contracts with the increase in value P"*". We note that it is 
not useful to use the spatial wave function to compute the form factor because of the appearance of the factor I — x 
in the exponential of Eq. 



VI. THE REST FRAME CHARGE DISTRIBUTION IS GENERALLY NOT OBSERVABLE 

The concept of a charge density that depends on three spatial variables, but not on the time, is inherently non- 
relativistic. This is because the use of only three variables involves replacing a four-dimcnsional quantity by one 
involving only three dimensions. One procedure, discussed above in Sect. II, is to evaluate the Feynman diagram of 
Fig. [l] by integrating over the component of the virtual momentum k. This leads to a formalism in which the 
form factor depends on a Fourier transform of the square of a wave function that depends both on position B and 
momentum x variables. 

One can try to recover the more familiar three-spatial dimension formalism by evaluating the Feynman diagram of 
Fig. [T] using the time-ordered-perturbation theory TOPT formalism in the rest frame. One proceeds by integrating 
over all times, with the exponential oscillating factors converted into energy denominators. In the TOPT formalism 
any given Feynman diagram is the sum of several TOPT diagrams. In the present case, the sum of the two TOPT 
diagrams of Fig. [5] leads to the Feynman diagram of Fig. [TJ Only Fig. [^h- corresponds to measuring a density. The 
term of Fig. [^h corresponds to the hadronic part of the incident photon wave function interacting with the target. 

One can examine the contribution of the term of Fig. [5^, F2a(Q^)- It is given Q by 

an expression that leads to the correct result in the infinite momentum frame (P — > oo) 0]. This expression can 
be interpreted as involving an initial and a final state wave function if one interprets the energy denominators 
(multiplied by phase space factors) as wave function expressed in momentum space. The symbol — > used above refers 



cc 




FIG. 3: Non-conservation of current is measured by the deviation of CC from as a function of = (in units of M^). 
Solid curve m = 0.51 B/M = 0.04, short-dashed curve m = 0.501 B/M = 0.004 long-dashed curve mO.6 B/M = 0.33 ( min 
units of M) 



m2 



to the approximation of keeping only a single TOPT diagram. For simplicity we take the example, mi 
and also work in the target rest frame: P = to isolate the rest frame charge distribution. Then Ep = V-P^ + = 
M, E-p+^ = ^(P + q)2 = v^q2 ^ M^, Ei = y^p^ + m^, E'l = ^(pTqpT^. The four vector (p^' +p"^) = 

[El -I- £"i, 2p -I- q]. There are three integrals appearing on the right-hand-side of Eq. (|48l) : 



qJ2(q') 
qJ3(q') ^ J 



(vV 



V(p + q)^ + "^^) 



2Ei2Ei'2E2 {Ep - Ei - E2){Ep+^ ~E\-E2) 
(fp 2p 
2Ei2Ei'2E2 {Ep - Ei - E2){Ep+^ -E\-E2) 

(fip q 
2Ei2Ei'2E2 {Ep - Ei - E2){Ep+q - E'l - E2) ' 



Thus we arrive at the four-vector equality 

{2P + qYF2a{Q'') 



{2^y 



r[/i,q(J2 + J3)], 



(49) 
(50) 
(51) 

(52) 



Maintaining current conservation requires that the matrix element of q^J^ vanishes. Taking the scalar product of 
Eq. ([5^ with leads to the requirement: 



- <7"/i - |q|(J2 + ,h) = CC /i(q2 = 0), 



(53) 



where go = -\/q2 + — M. The right-hand-side of Eq. ((53)) is defined as CC /i(0) so that comparing CC to unity 
provides a reasonable measure of the failure of this approximation to uphold current conservation. We express all 
momenta and mass in units of the target mass (=1), take as an examples m — 0.501, 0.51 and 0.6 and plot the 
numerical results in Fig. [31 Both CC and q^ — g^j-g measured in units of the target mass M, taken as unity. Thus 
the natural scale of any quantity is unity. We see that current conservation is massively violated in the rest frame for 
systems in which B/m = (2m ~ M)/ m is not very small. In that case, the expression that potentially depends on the 
square of the wave function or density has no independent physical reality. 
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FIG. 4: Form factor for a neutral system with one heavy m2 = M and one light mi = 0.14 M negatively charged constituent. 



VII. NEUTRAL SYSTEMS 

Previous work [sl showed that the central transverse density of the neutron is negative. This contrasts with the 
long held view that there must be positive charge density at the center to neutralize the effects of a negatively charged 
pionic cloud that occupies the exterior. This result demands interpretation [l6l|-|l8|. 

One relevant question is whether or not the intuition that a neutral system consisting of a heavy charge positively 
charged particle and a negatively charged lighter particle disobeys the standard intuition that the averaged squared 
charged radius is negative, when the charge density is evaluated in the infinite momentum frame. We examine this 
question in our model by taking the (of mass mi) to be positively charged and the ^ (of mass m2 < rrii to be 
negatively charged. 

The form factor of this model can be obtained by using Eq. ([8]) by including a second term obtained by interchanging 
mi and TO2 and putting a minus sign in front. That operation gives the result 



F{Q') = ^ / dxx 



Tanh-i^ VQ^(i-) 



47r^ Jo \fW\J^x ra\ + 4(1 - x)m\ - x{\ - x)M'^ + (1 - xfQ^ 

^^42: mf+4(l-i:)m|-a;(l-i:)Af2 + (l-a;)2Q2 J ^ 



^jQ^^J^x m\ + 4(1 - x)'m\ - x(l - x)^^ + (1 - xfQ^ ' 

The results of a numerical evaluation using mi — M and m2 — 0.14M are shown in Fig. [4l One observes the rise of 
F{Q'^) from zero, which is the effect expected from non-relativistic, rest-frame considerations. The effective squared 
radius, defined in Eq. ()43|) is indeed negative. This is the same as expected from the intuition that the negatively 
charged light particle resides on the outside edge of the system. 

One obtains the analytic result for the charge radius by taking the limit of very low in the expression for the 
form factor Eq. (|54p . The expression is simplified if one uses the (relevant for nucleon)[l9| case of mi = AI. Then 



11 





bM 



FIG. 5: Transverse charge density for a neutral system of a positively charged heavy object and a negatively charged lighter 
object 



one finds: 




The radius is dominated by a singular term proportional to 1 /m\ . Thus as expected the lighter constituent drifts to 
the edge of the nucleon. The conventional expectation is borne out on the light front. This is shown in more detail by 
plotting bp{b), as shown in Fig. [5] The positive charge density is concentrated at the center and the negative at the 
edge. This finding does not contradict the explanations offered in Refs. [l^-[l5. Ref . fl6l argues that negative charge 
at high X corresponds to negative charge at small values of b. The Nn model of Ref. shows that one must include 
the finite size of the nucleon to obtain a computed Fi that looks like the measured function. Thus the point-like 
nature of the constituents used here is unsurprisingly not relatistic. Moreover, in that model negatively charged pions 
reside both at the edge and at the center of the nucleon. The implication of [l^ is that the pions may have large 
values of longitudinal momentum fraction. This expectation is borne out by the model calculation [20| . Thus in pion 
cloud models of the nucleon pions that have a large longitudinal momentum tend to reside near the center of the 
nucleon. 



VIII. NON-RELATIVISTIC LIMIT 



The conventional lore is that the electromagnetic form factor is the Fourier transform of the charge density. In this 
section we see how this idea emerges by taking the non-relativistic limit. 

Our starting point is the wave function Eq. (fT6| and the form factor Eq. p7)) . Recall that the quantity x = k'^ / . 
We work in the rest frame and take the non-relativistic limit in which the energy /c'^ = mi, and fc+ = mi + , where 
is the third-component of the relative longitudinal momentum. Then [2^ . [25| 

mi -\- M — mi — TO2 — B — , . 

X — — — , l — x— — — , (57) 

where in conformation with non-relativistic notation, we define the positive binding energy B so that 



M = mi + 7712 — B. 



(58) 
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+ "^2 




X 


1 


— X 




= 


- M\ 


K,^ + m\ 


+ mn , 




mi + 


m2 — B — 






- M\ 


K? + m\ 


1 + — f 

\ mi mi / 


+ mi 
+ 


mi 


TO2 



To obtain the non-relativistic expression we express the denominator appearing in Eq. (|16p in terms of k^. This gives 

(59) 
(60) 

+ + ' 61 

m2 TO2 / 

^ {n^ + kI){— + —) + mi+m2 + B] (62) 
mi m2 

= 2M(-i?-^), (63) 

where 

9 9 ■> ^ m^mo 
k'^ = + ni, K = K + K-^z, u = . (64) 

TOi + TO2 

In going from Eq. (l59|) to Eq. ([63|) we have ignored terms in v/c = k/m oi order three and higher. The result is that 
Eq. ((63|) is recognizable as 2M times the inverse of the non-relativistic propagator. 

The next step is to determine the coordinate form of the non-relativistic wave function ipNRi'r) (where r is canonically 
conjugate to k) and to show that the non-relativistic form factor is a three-dimensional Fourier transform of \iPnr{'P)\^ ■ 
First use the non-relativistic approximation Eq. ((63| in Eq. (fT6|) to find 



^NRi^i) = A' = 2fiB. (65) 

The coordinate-space wave function ^NR(f) is given by 



The expression Eq. (j66p is seen as the standard result obtained for the bound state of a two-particle system interacting 
via an attractive delta function potential. 

The wave functions Eq. (|65p and Eq. (|66p enable us to examine the condition needed for the approximations Eq. fST]) 
to be valid. For Eq. (|57|) to work we need <C mf 2, but from the wave functions ~ so that we require 

^ « 1 (67) 
"^1,2 

for the non-relativistic approximation to be valid. More specifically let M be the lighter of mi,m2, then wc may 
write the approximate condition as 

^ « 1- 

The non-relativistic form factor Fnr{Q^) is obtained by using Eq. (p5|) in the expression for the form factor Eq. (|17p . 
and taking the non-relativistic limit defined by the expressions: 

dx ^ (69) 

(mi + TO2) 

^(1 _ ^) _ (70) 
mi + m2 

(i-^)q^ 7 q- (71) 

mi -I- m2 



The result is 

2(27r)3/i 



Fnr{Q') = t^tAt- I d'r |^jvfl(r)|' e''^-"^^ . (72) 
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This is the usual expectation that the form factor is a three-dimensional Fourier transform of the wave function. We 
may evaluate the integral immediately to find 

771 /ri^\ 2(mi+m2)A /„„x 

2(mi+m2) A 

where Q = |q| and the coupling constants and other constants enter in such a manner as to make Fnr{Q'^ = 0) = 1. 

In the remainder of this section we study the accuracy of the non-relativistic approximation by comparing the 
results of using Eq. ((75|) with the model-exact results of using Eq. ^ for several examples. 



A. Bound state of two equal mass particles 



With equal masses mi = TO2 = ni the bound state can be thought of as a toy meson or a deuteron. Use mi m2 = m 

(2) / 
NR\ 



in Eq. (|65p leads to the non-relativistic wave function with 



with 



Xl = mB. (75) 

The coordinate space wave function tpNRir) is then 

Thus the wave function is the usual bound state wave function one obtains with a delta function binding interaction. 
We obtain the non-relativistic version of the form factor by using mi — m2 = m, A ^ A2 in Eq. (|73p to find 

tan~i 

FnUQ') = (77) 
4A2 

where Q = |q| and the coupling constants and other constants enter in such a manner as to make fI^I{Q^ = 0) = 1. 



We study the non-relativistic approximation numerically by comparing the exact model results Eq. (jl7[) with those 
of the non-relativistic approximation Eq. (j77p . See Fig. [Sj The figure shows two sets of results. In the upper panel 
the binding energy B = 0.002 M. This corresponds roughly to deuteron kinematics, in which the binding energy 
is of the order of a 0.004 of the deuteron mass. We see that the non-relativistic approximation is not accurate for 
values of Q^/M^ greater than about 1. If one increases the binding energy to 0.1 AI, one sees that the non-relativistic 
approximation is not accurate for any value of Q^. If one approximates a nucleon by taking Al — 1 GeV, then 
m = 0.55 GeV, which is much larger than a, u,d constituent quark mass. Thus the range of masses for which the 
non-relativistic approximation is valid is very narrow indeed. 

We can gain some insight into the nature of the relativistic corrections to the charge radius by studying the low 
limit of the form factor of Eq. ([5]). One finds 

/n2 D*2 

where we use the notation R*^ to denote an effective radius squared that is not generally associated with the expec- 
tation of the square of a radius operator weighted by a density. The explicit evaluation gives 

, , (4+487)cot-i(27)+ 2 _24 

M^R*^ = ^^-^^ ^ ^ ^, (79) 

16((27+2^)cot-i(27)-l 



and 
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FIG. 6: Exact vs non-relativistic form factors for the case mi = m2 = m. Solid curve- exact result, dashed curve-non-relativistic 
limit. Upper panel deuterium-like kinematics in which B = 0.002Af. Lower panel- B = O.IM. 



The non-relativistic limit corresponds to the limit of small values of 7, which corresponds to a small value of B/M . 
So we expand the previous result to order B/M to find 



„ ^2 (12288 - 28167r2 -f 195^4) 5 (128^2 - 25\/27r2) 64 - 5^^ ^2 Af 
^^^^^ 48M^^ i^r^ + + + ^''^ 



The non-relativistic value of the mean square radius, Rjfji (which is a true mean square radius) is obtained by 
expanding the form factor for small values of Q^: 

" aWb' ^^^^ 



which corresponds to the leading term of Eq. (|5T|) in the limit that B approaches 0. Comparing Eq. (|5T|) with Eq. 
shows that the former contains a series of terms that represent the boost corrections to the non-relativistic result. 
Each correction is positive and can be substantial. The figure [7] shows the ratio of the exact of the mean square radius 
to the non-relativistic approximation as a function of B/M. We see that the non-relativistic approximation works 
well only for very small values of B/M. Indeed, the ratio of the leading correction to the non-relativistic result is 
given by 

For this ratio to be less than 10 %, must be less than one part in a thousand! Thus, within the framework of our 
toy model, the relativistic corrections can generally expected to be very substantial. 



B. Quark-diquark model of the nucleon 

Another interesting example is motivated by recent quark-diquark models of the nucleon [2l|, [l^, H^. We take 
mi = m, 1712 = 2mi = 2m. Then from Eq. (|58p we have M = 3m — B. In these, models current quarks acquire a large 
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FIG. 7: Ratio of exact to non-relativistic effective square radii for the case mi — 7712 = m as a function of the ratio of the 
binding energy B to the hadronic mass M. This is also the ratio of the true value of (6^) to its non-relativisitic version. 



constituent mass due to the effects of dynamical chiral symmetry breaking. Therefore we take m = 400 MeV, and 
M = 940 MeV which corresponds via Eq. ((551) to B = 260 MeV and B/M = 0.276. The non-relativistic expression 
for the form factor, Fff^(Q^) for this case is obtain using the appropriate reduced mass as 



n'^m - (84) 

3Ai2 

A?2 = t'^B. (85) 



3 

Results comparing the exact form factor computed from Eq. ([8]) with that of Eq. ((84)) are shown in Fig. [8l The non- 
relativistic version gives a poor approximation to the exact form factor for all values of Q^. This can be understood 
by considering the effective squared radius R*q2q for this case. We find 



I£j2q^R%R 6 -In 2 

J L jv 



(86) 



"-NR. 



The right-hand-side is evaluated as 0.887 for the present case, so that there is a substantial relativistic correction to 
the quantity i^j^jj — 1 for any non-zero value of Q^. This means that one can not take a three-dimensional Fourier 
transform of the form factor to get a charge density even if the constituent masses are large. 



C. Nuclear Physics and mi 7^ m2 



We consider masses that correspond to electron scattering from a charged nucleon of mass m (which is the free 
nucleon mass minus the average binding energy per nucleon of 8 MeV) bound in a nucleus of mass M = mA, with a 
spectator system of mass m2 = {A — l)m -I- S, where S is the orbital separation energy. We measure all momenta in 
terms of m = 932 MeV, and take the separation energy S — 0.05 or Sm about 46 MeV. The results for A = A and 
A = 208 are shown in Figs. [9] and [TOl The startling finding is that the relativistic effects reduce the form factor for 
light nuclei, but increase it for heavy nuclei. Furthermore, the relativistic effects are larger for heavy nuclei than for 
light nuclei (for a fixed value of S.) 

We obtain some analytic understanding by expanding the effective squared radius (defined in Eq. (j78p in powers of 
S. We find 



^ ^ vglTP (4A - 2) log ((^ - 1)')) 
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FIG. 8: Exact vs non-relativistic form factors for the case m2 = 2mi,m = 400 MeV, B = 260 MeV = 0.276 M. Solid 
curve-exact, dashed non-relativistic. 




FIG. 9: Exact vs (solid curve) non-relativistic form factors (dashed curve) for A — 4 for Sm = 46 MeV. 

We see that the first term is indeed the non-relativistic result, and that the second term changes sign for the value of 
A that satisfies the equation iA — 2{A — 2) ln{A — 1) = or A w 12. This displayed in Fig. [TT] It is also seen that the 
relativity causes very significant effects on the effective radii. Except for values of A near 12, the changes are of the 
order of 10-15%. I expect that the specific values shown in Fig. [Tl] are highly model-dependent. Covariant models 
other than the 5*0^ model used here probably have have effects of different sizes. However, the large effects shown 
here cause one to wonder if relativity really may cause the true nuclear radii extracted from elastic electron scattering 
to differ by 10-20% from those appearing in tables. As noted above, we can expect that the model employed here is a 
reasonable representation of the lowest s-state of heavy nuclei for which the range of the binding interactions is much 
less than the size of the system as a whole. For such states, the results of Fig. [TT] should be a reasonably accurate 
guide, so that significant effects of relativity should be expected. 
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FIG. 10: Exact (solid curve) vs (dashed curve) non-relativistic form factors for A = 208 for Sm = 46 MeV. 




FIG. 11: Exact vs non-relativistic effective radii, -M — j- — 1 as a function of A for Sm — 46 MeV. This is also the ratio of {b"^) 
to its non-relativistic counterpart. 

IX. SUMMARY 

A relativistic model of a scalar particle 4" is a bound state of two scalar particles </> and ^ is used to elucidate 
relativistic aspects of electromagnetic form factors. First, the form factor for the situation in which the ^ and 4> 
carry a single unit of charge, but the ^ is neutral is computed using an exact covariant calculation of the lowest-order 
triangle diagram. This is followed by a another derivation using the light-front technique of integrating over the 
minus-component of the virtual momentum in Sect. Ill that obtains the same form factor. This is also the result 
obtained originally by Q by using time-ordered perturbation theory in the infinite-momentum-frame IMF. Thus three 
different approaches yield the same exact result for this model problem. The asymptotic limit of asymptotically high 
momentum transfer is also studied with the result that F{Q^) ~ 1/2 In^ /Q^ The next section (IV) explains the 
meaning of transverse density p{h) of the model. Its central value varies singularly as ln'^(6)/3. A general derivation 
of the relationship of p{h) with the form factor using three dimensional spatial coordinates is presented. This allows 
us to identify a mean-square transverse size (6^) = J (Pb b^p{b) that is given by &^ = —^-^t{Q^ = 0). The quantity 
(6^) is a true measure of hadronic size because of its direct relationship with the transverse density. Using this model 
it is possible to display the spatial wave function in terms of three spatial coordinates (Section V), but this is not 
very useful. Section VI shows that the rest-frame charge distribution is generally not observable by studying the 
explicit failure to uphold current conservation. Section VII shows that neutral systems of two constituents obey the 
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conventional lore that the heavier one is generally closer to the transverse origin than the lighter one. It is also 
argued that the negative central charge density of the neutron arises in pion-cloud models from pions residing at 
the center of the nucleon. The non-relativistic limit is defined and applied to a variety of examples in Section VIII. 
By varying the masses one can study a continuum of examples in which the constituents move at a wide range of 
average velocities. The relevant quantity is the ratio of the binding energy B to that of the mass Ai of the lightest 
constituent (0 or <^). For small values oi B/M the exact relativistic formula is shown to be the same as the familiar 
one of the three-dimensional Fourier transform of a square of a wave function. If the and ^ have equal masses m 
we find that B/{2m) must be less than 0.001 for the relativistic corrections to mean-square radii to be be less than 
10%, see Eq. ((83)) . For the case when = 2m^ which mimics the quark-di-quark model of the nucleon we find that 
there are substantial relativistic corrections to the form factor for any value of Q^. This means that one can not take 
a three-dimensional Fourier transform of the form factor to get a charge density even if the constituent masses are 
large. A schematic model of the lowest s-states of nuclei is developed by choosing = {A — l)m^, where A is the 
nucleon number. Relativistic effects are found to decrease the form factor for light nuclei but to increase the form 
factor for heavy nuclei. Furthermore, these lowest s-states are likely to be strongly influenced by relativistic effects 
that are order 15-20%. 

I thank the USDOE (FG02-97ER41014) for partial support of this work, S. Brodsky for advocating the use of the 
0"^ model as a pedagogic tool, and J. Arrington, A. Bernstein, M. Burkardt, I. Cloet, B. Jennings, E. Henley, and W. 
Polyzou for useful discussions. 
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